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I. INTRODUCTION 



Lepton number (L» = L e , L^, L T ) is conserved in the standard model (SM) with massless 
neutrinos. However experimental data such as those coming from Super-Kamiokande [l[ 
Q or Sudbury Neutrino Observatory 3] showed evidence of neutrino oscillations. These 
experimental results imply lepton flavor violation (LFV) and the lepton sector of the SM 
requires an extension including massive neutrinos. Since conservation of the lepton flavor 
number is not a law of nature that we rely on, it is significant to quest for new physics via 
lepton flavor violating effective operators. Effective lagrangian extension of the SM could 
have an impact on particle physics and astrophysics. Therefore it is crucial to investigate 
the physics potential of neutrino experiments to probe LFV. 

Neutral-current lepton flavor violating processes such as Z — > ifij = e, /x, r) can 
be observable at colliders and effective operators contributing to these processes have been 

nn 

stringently constrained by collider experiments |4J, |5| . On the other hand, in colliders neu- 
trinos are not detected directly in the detectors. Instead, their presence is inferred from 
missing energy signal. Therefore it is impossible to detect a neutrino flavor at a collider 
and charged- current processes such as W ± — > ifuj can not be discerned. In order to con- 
strain lepton flavor violating couplings of neutrinos to W boson, beta beams provides an 
excellent opportunity. Beta beams are electron neutrino and antineutrino beams produced 
via the beta decay of boosted radioactive ions [6j. Such decays produce pure, intense and 
collimated neutrino or antineutrino beams. Purity of the neutrino beam make it possible to 
probe couplings WiiV e with a high precision. 

In the original beta beam scenario ion beams are accelerated in the proton synchrotron or 
super proton synchrotron at CERN up to a Lorentz gamma factor of 7 ~ 100, and then they 
are allowed to decay in the straight section of a storage ring. After the original proposal, 
different options for beta beams were investigated. A low gamma (7 = 5 — 14) option was 
first proposed by Volpe [7]. Physics potential of low-energy beta beams was discussed in 
detail. It was shown that such beams could have an important impact on nuclear physics, 
particle physics and astrophysics Higher gamma options for the beta beams have 

also been studied in the literature 14J, |20|429|]. A higher gamma factor provides several 
advantages. Firstly, neutrino fluxes increase quadratically with the gamma factor. Secondly, 
neutrino scattering cross sections grow with the energy and hence considerable enhancement 
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is expected in the statistics. An additional advantage of a higher gamma option is that it 
provides us the opportunity to study deep inelastic neutrino scattering from the nucleus. 
Very high gamma (~ 2000) options would require modifications in the original plan such 
as using LHC and therefore extensive feasibility study is needed. In this context medium 
energy setup is more appealing and less speculative. 

In this paper we investigate the physics potential of a medium energy setup (7 = 350 — 
580) proposed in Ref. [20j to probe lepton flavor violating couplings W \iv e . As far as we 
know, phenomenology of LFV in a beta beam facility was studied only in Ref. [30] but 
considering supersymmetric models. Different from [30] we have probed LFV in a model 
independent way by means of the effective Lagrangian approach. 

The organization of this paper is as follows: In the next section we outline the effective 
Lagrangian approach. In section III we summarize the neutrino fluxes and the cross sections 
for quasielastic and deep inelastic scattering and present our main results. Finally Section 
IV includes concluding remarks. 



II. EFFECTIVE LAGRANGIAN FOR LEPTON FLAVOR VIOLATING WliVj 
COUPLINGS 



There is an extensive literature on non-standard interactions of neutrinos |29|, |31|-|40]. 
New physics contributions to lepton flavor violating WliVj couplings can be investigated in 
a model independent fashion by means of the effective Lagrangian approach. The theoretical 
basis of such an approach rely on the assumption that at higher energies beyond where the 
SM is valid, there is a more fundamental theory which reduces to the SM at lower energies. 
The SM is assumed to be an effective low-energy theory in which heavy fields have been 
integrated out. Such a procedure is quite general and independent of the model at the new 
physics energy scale. 

Specifically we consider the SU(2)l <8> U(l)y invariant effective Lagrangian introduced in 
Refs. Q, 41]. The effective Lagrangian can be written as 

C eff = C SM + i + ... (dim > 6) (1) 

n,i,j 

where i,j = 1, 2, 3 denote flavor indices, n runs over the number of independent operators of 
dimension-6 and A is the energy scale of new physics. c%[ are the dimensionless anomalous 



coupling constants. There are four independent operators of dimension-6 contributing to 
WiiVj vertex: 



O l l w = CUl»r a D»L J )W* v + h.c. (2) 
0f L 13 = i{tpT*D ll tl>){L(fT»L j ) + h.c. (3) 
0% e = (L i Dp£ Rj )D»<f> + h.c. (4) 

0%^ = L i( j^T a £ Rj <f>W^ + h.c. (5) 

where and £rj are the left-handed lepton doublet and right-handed singlet of SU(2)l <S> 
U(1)y- is the scalar doublet and is the covariant derivative. W% v and r a are the 
SU (2) L gauge boson field tensors and Pauli matrices respectively. 

After spontaneous symmetry breaking, operators ([2H5D give rise to vertex functions for 
WiiVj. The vertex functions for W (k)£i(p2)uj(pl) generated from the effective Lagrangian 
are given, respectively, by 

rf 1} = Cl V2 [-p x k^P L + P$k^P L - (k ■ P2 )i X Pl + (k ■ Pi)l x P L } (6) 

rf 2) = ~ c ^ Xp L (7) 
r? 3) = -<*f v\Pl (8) 

T x {4) = 2ir 1 c A k p a px P L (9) 



For a convention, we assume that p\ is incoming to and p 2 and k are outgoing from the 

c (3) 

vertex, q (i = 1, ..,4) are scaled coupling constants defined by: c\ = £ j^f, c 2 = -jt, C3 = 
and C4 = (For abbreviation we drop flavor indices). Pl = |(1 — 75), g is SU(2)l 

coupling constant and 77 represents the vacuum expectation value of the scalar field. (For 
definiteness, we take 77 = 246 GeV in the calculations presented in this paper). 

Operators (J2j(5]) not only contribute to WiiVj but also Z£~£^ couplings. On the other 



hand, Z£ { £~j~ receive contributions from 9 independent operators [41|. Therefore processes 
involving Z£~£^ vertex do not isolate the contributions coming from operators ([2}|5]). 
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III. NEUTRINO FLUXES AND CROSS SECTIONS 



In a beta beam facility very intense and collimated neutrino or antineutrino beams can 
be produced by accelerating /3-unstable heavy ions to a given 7 factor and allowing them to 
decay in the straight section of a storage ring. In the ion rest frame the neutrino spectrum 
is given by the following formula 

dN 



dcosOdE, EHK-EMEr-W-nt (10) 

where Eq is the electron end-point energy, m e is the electron mass. E u and 9 are the energy 
and polar angle of the neutrino. Neutrino flux observed in the laboratory frame can be 
obtained by performing a Lorentz boost. The neutrino flux per solid angle in a detector 
located at a distance L is given by [2o| 

where < y < 1 - y e , y = y e = and 

Ve 




9(Ve) = ™ VI - l/e 2 (2 - 9y 2 e - 8y 4 e ) + 15ytLog 



(12) 



18 Ne and 6 He ions have been proposed as ideal candidates for a neutrino and an antineutrino 
source, respectively 0, Q]. These ions produce pure (anti)neutrino beams via the reactions 
l^Ne — >l s Fe + v e and \He ++ — y\ Li +++ e~u e . We assume that total number of ion decays 
per year is Np = 1.1 x 10 18 for 18 A^e and N p = 2.9 x 10 18 for 6 He. 

Neutrino and antineutrino fluxes as a function of (anti) neutrino energy at a detector of 
L = 732 km distance are plotted in Fig. [TJ 7 parameters for ions are taken to be 7 = 350 
for 6 He and 7 = 580 for 18 iVe. The foregoing detector distance and 7 values have been 



proposed in Ref. |20j as a medium energy setup. In Ref. |20j authors have considered a 
Megaton-class water Cerenkov detector with a fiducial mass of 400 kiloton. They show that 
a cut demanding the reconstructed energy to be larger than 500 MeV suppresses most of 
the residual backgrounds. We assumed a water Cerenkov detector with the same mass and 
a cut of 500 MeV for the calculations presented in this paper. 

In Fig{T]we observe that neutrino spectrum extend up to 4 GeV and antineutrino spectrum 
extend up to 2.5 GeV. The energy range of the neutrino spectrum is comparably larger than 
the antineutrino spectrum. Therefore number of events for antineutrinos is expected to be 



low. So we do not perform an analysis for antineutrinos. For neutrino energies between 0.5 
-1.5 GeV, dominant contribution to the cross section is provided by quasielastic scattering. 
When neutrino energy exceeds 1.5 GeV, deep inelastic scattering starts to dominate the cross 
section. Neutrino electron scattering takes place at all energies in the spectrum. Lepton 
flavor violating v e e~ — > v e yT process contains both W and Z exchange diagrams. Hence 
this process receives contributions from both W fiv e and Z[ie couplings. Therefore it is not 
possible to set limits on W ' [iv e coupling independent from Z[ie. Since the main advantage 
of a beta beam facility has been lost (isolation of the W /iz/ e vertex) we do not analyse the 
process v e e~ — > v e \i~ . 



A. Neutrino nucleon quasielastic scattering 



In addition to SM quasielastic scattering v e n — > e~ p, operators flH|5]) give rise to new 
reactions v e n p and v e n — > r~ p. The tau production threshold is 3.5 GeV. Hence, 

tau production calls for a higher 7 factor and medium energy setup of a beta-beam facility 
is not convenient to analyse v e n — > r~p reaction. On the other hand, muon production via 
quasielastic scattering seems to be appealing. Feynman diagram for v e n — > p is given 



in the left panel of 



standard formula 42| . 



'igfSJ We see from the diagram that hadron current is given by the 
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2V2 



cos 6c u p 



1 



2m N 



w 



U r 



(13) 



where cosOc = 0.974 is the Cabibbo angle, is the mass of the nucleon and F's are 
invariant form factors that depend on the transferred momentum q 2 = (p p — p n ) 2 . The F's 
are known as vector Fy, axial- vector Fa and tensor F\y (or weak magnetism) form factors. 
They are all G-parity invariant. We adopt the same parameterization of the momentum 
dependence as in Ref. [431 ] : 



Fy(q 2 ) 
Fw(q 2 ) 
F A (q 2 ) 



4ml 



1 - 



4m, 



(0.84 GeV) 2 



4mi 



(14) 



1.270 1 - 



(0.84 GeV)'\ 
2 \ -2 



Q 



1.032 GeV) 2 



Here £ = (/i p — ji n ) / Hn = 3.706 is the difference in the anomalous magnetic moments of the 
nucleons. Lepton current generated from the effective vertices dMH]) is given by 



j* = u, [v x {1) + r* , + rj 3) + rf 4) ] u Ve (15) 

From ( !T3|) and ( fl5]) scattering amplitude is obtained as 



M = -S sP" ( 16 ) 

q 2 -m 2 w 

where mw is the W boson mass. The W propagator can be approximated as; q 2_™z v ~ 

In FigfS] we show total cross section of v e n — > fi~ p as a function of neutrino energy 
for some values of the anomalous couplings. We see from the figure that cross sections 
proportional to couplings c\, C3 and C4 have similar behaviors as a function of neutrino 
energy. They all increase as the energy increases but increment rate is especially high 
up to 2 GeV. On the other hand, anomalous cross section for C2 attains its maximum at 
approximately 1 GeV and then starts to decrease. 



B. Charged-current deep inelastic scattering 



When neutrino energy exceeds 1.5 GeV, deep inelastic scattering starts to dominate 
the cross section. Charged-current deep inelastic scattering of an electron-neutrino from 
nucleon is described by t-channel W exchange diagram (right panel of FigfS]). Since quark 
couplings to W boson are not modified by operators (EED hadron tensor does not receive 



any contribution. It is defined in the standard form 



44 1 



W, 



q z J p ■ q 



fii/a/3 



qOp/3 

2p ■ q 



F 3 (x,Q 2 



(17) 



where p^ is the nucleon momentum, q^ is the momentum of the W boson propagator, 



Q 2 = -q 2 , x 



2p-q 



and 



. _ P'Q 



The structure functions for scattering on a proton are defined as follows [44 1 



2x(d + u + c + s) 
2(d — u — c + s) 



(18) 



The form factors Fi's can be obtained from (1181) by using Callan-Gross relation 2xF\ = F 2 



45(. The structure functions for scattering on a neutron are obtained from those of the 
proton by the interchange u •<-)■ d. In our calculations parton distribution functions of 
Martin et a/.|46j have been used. We assumed an isoscalar oxygen nucleus iV = (p + n)/2 
and two free protons for each H 2 molecule. Naturally occurring oxygen is 99.8% 16 which 



is isoscalar |47|. Hence the error incurred by assuming an isoscalar oxygen target would be 
not more than a fraction of one percent. 

Possible new physics contributions coming from the operators (J2H5J) only modify the lepton 
tensor: 

V = Hk')T,u{k)} ] [u{k')V v u{k)} (19) 

spin 

where k and k! are the momenta of initial v e and final fi~, respectively. T u represents the 
sum of anomalous vertices (jMH])- 

The muon production via charged- current deep inelastic scattering of neutrinos from the 
proton is plotted in FigJH We see from the figure that different from quasielastic scattering 
case, increment rate is higher at high energies. 

C. Statistical analysis and results 

A detailed investigation of the anomalous couplings requires a statistical analysis. To this 
purpose, number of events have to be calculated. If we assume that initial neutrino beam is 
pure and contains only electron-neutrinos then SM cross section for muon production is zero. 
But neutrino beam near the detector should contain a small fraction of muon-neutrinos due 
to neutrino oscillations. In Figj5] we present the transition probability P[y e — > z/ M ) and the 
survival probability P[y e — > v e ) as a function of neutrino energy at a detector of distance 



732 km. We use the following approximate formulas 48] 



Piye \ sm 2 (2fl 13 ) sin 2 6 23 (l - cos (^2.54^^ ) ) (20) 

P(i/. ->!/.)= l-^sin 2 (2M(l- cos (2.54^^)) (21) 

where Am| 3 is the neutrino mass-squared difference in eV 2 , L is the distance between neu- 
trino source and detector in m and neutrino energy E v is given in units of MeV. In the 
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calculations we assume that sin 6 13 = 0.035, sin 2 023 = 0.50 and |Am|,| = 2.40 x 10~ 3 eV 2 



491 ] . Number of events has been obtained by integrating cross section over the neutrino 



energy spectrum and transition probability and multiplying by the appropriate factor that 
accounts for the number of corresponding particles (protons or neutrons) in a 400 kiloton 
fiducial mass of the detector. For instance, number of SM events for charged-current deep 
inelastic scattering is given through the formula, 

N SM = f P{Pe p„) (Jgj]f) [ N p^p^-x + N n <T Vll «r+p-x]dE v (22) 

where Uu^p-^^-x and cr v n -*n-x are the SM cross sections of deep inelastic scattering of the 
muon-neutrino from the proton and neutron respectively. N p and N n are the number of 
protons and neutrons in a 400 kiloton fiducial mass of the detector. Number of lepton flavor 
violating events can be calculated in a similar manner but considering lepton flavor violating 
cross sections and survival probability. 

We studied 95% confidence level (C.L.) bounds using one-parameter \ 2 analysis without 
a systematic error. The y 2 function is given by, 

-' = {n\-]' < 23 > 

\1\SM Ostat J 

where N is the number of lepton flavor violating events as a function of couplings ci, €2,03,04, 
Nsm is the number of events expected in the SM and S sta t is the statistical error. In Table 
[U we show 95% C.L. bounds on the couplings Ci, 02,03 and C4 obtained from quasielastic 
and deep inelastic scatterings. These bounds are obtained for 1 year running time of the 
beta-beam experiment. We see from the table that bounds obtained from quasielastic and 
deep inelastic scatterings are close to each other. Therefore, they have almost same potential 
to probe LFV. 



IV. CONCLUSIONS 



Beta beams present an ideal venue to measure neutrino cross sections. For beta beams 
neutrino fluxes are precisely known and therefore uncertainties associated with the neutrino 
(antineutrino) fluxes are negligible. Purity of the produced (anti) electron-neutrino beam is 
an other advantage of beta beams. These features enable to detect neutrino cross sections 
with a high precision. 
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Medium energy setup of a beta-beam facility provides us the opportunity to isolate W/j,u e 
vertex which is not the case for a collider experiment. Probing W \w e couplings is important 
for understanding the physics beyond the SM and contributes to the studies in neutrino 
physics. In neutrino oscillation experiments, identification of the neutrino flavor is based 
on charged-current processes at the detector. Neutrino flavor can be identified through 
the flavor of the associated charged lepton. Hence, lepton flavor violating couplings of the 
neutrinos to W can mimic the oscillation signals. Therefore constraints on W \w e couplings 
are important for precision measurements of the oscillation parameters. 
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FIG. 1: Fluxes as a function of neutrino energy for v e (solid line) and v e (dotted line) at a detector 
of L = 732 km distance. 7 parameter is taken to be 350 for v e and 580 for v e . 



w 




FIG. 2: Figure on the left shows Feynman diagram for u e n — > fi p. Figure on the right represents 
schematic diagram for neutrino deep inelastic scattering u e N — > yTX. 



TABLE I: 95% C.L. bounds on the couplings 01,02,03 and C4 obtained from quasielastic and deep 
inelastic scatterings. Only one of the couplings is assumed to be non-zero at a time. 



Quasielastic scattering: 










|ci| < 0.0159 


|c 2 | < 8.457 x 10~ 7 


|c 3 | < 0.0010 |c 4 | < 8.561 x 10~ 5 


Deep inelastic scattering: 










|d | < 0.0105 


|c 2 | < 8.902 x 10~ 7 


|c 3 | < 0.0010 |c 4 | < 7.846 x 10" 5 
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FIG. 3: Total cross section of v e n — >■ fi p as a function of neutrino energy. The values of the 
anomalous couplings is stated in the figure. 
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FIG. 4: Total cross section of charged-current deep inelastic scattering u e p — > n X as a function 
of neutrino energy. The values of the anomalous couplings is stated in the figure. 
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FIG. 5: Figure on the left is the transition probability P(y e — > v^) and figure on the right is the 
survival probability P(y e — > v e ) in vacum as a function of neutrino energy. Detector distance is 
taken to be 732 km. 
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